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TECHNIQUE FOR BIT-ACCURATE FILM GRAIN SIMULATION 

CROSS-REFERENCE TO RELATED APPLICATIONS 

5 This application claims priority under 35 U.S.C. 1 19(e) to U.S. Provisional Patent 

Application Serial No. 60/511,026, Filed on October 14, 2003, the teachings of which are 
incorporated herein. 

TECHNICAL FIELD 

10 

This invention relates to a technique for simulating film grain in an image. 
BACKGROUND ART 

15 Motion picture films comprise silver-halide crystals dispersed in an emulsion, which 

is coated in thin layers on a film base. The exposure and development of these crystals form 
the photographic image consisting of discrete tiny particles of silver. In color negatives, tiny 
blobs of dye occur on the sites where the silver crystals form following chemical removal of 
the silver during development of the film stock. These small specks of dye commonly bear 

20 the label •grain' in color film. Grain appears randomly distributed on the resulting image 
because of the random formation of silver crystals on the original emulsion. Within a 
uniformly exposed area, some crystals develop after exposure while others do not. 

Grain varies in size and shape. The faster the film, the larger the clumps of silver 
formed and blobs of dye generated, and the more they tend to group together in random 

25 patterns. The term "granularity" typically refers to the grain pattern. The naked eye cannot 
. distinguish individual grains, which vary from 0.0002 mm to about 0.002 mm. Instead, the 
eye resolves groups of grains, referred to as blobs. A viewer identifies these groups of blobs 
as film grain. As the image resolution becomes larger, the perception of the film grain 
becomes higher. Film grain becomes clearly noticeable on cinema and High Definition (HD) 

30 images, whereas film grain progressively loses importance in Standard Definition (SD) and 
becomes imperceptible in smaller formats. 

Motion picture film typically contains image-dependent noise resulting either from 
the physical process of exposure and development of the photographic film or from the 
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subsequent editing of the images. Photographic film possesses a characteristic quasi-random 
pattern, or texture, resulting from physical granularity of the photographic emulsion. 
Alternatively, simulation of similar pattern can occur in computed-generated images in order 
to blend them with photographic film. In both cases, this image-dependent noise bears the 
5 designation of "film grain." Quite often, moderate grain texture presents a desirable feature 
in motion pictures. In some instances, the film grain provides visual cues that facilitate the 
correct perception of two-dimensional pictures. Film grade often varies within a single film 
to provide various clues as to time reference, point of view, etc. Many other technical and 
artistic demands exist for controlling grain texture in the motion picture industry. Therefore, 

10 preserving the grainy appearance of images throughout image processing and delivery chain 
has become a requirement in the motion picture industry. 

Several commercially available products have the capability of simulating film grain, 
often for blending a computer-generated object into natural scene. Cineon® from Eastman 
Kodak Co, Rochester New York, one of the first digital film applications to implement grain 

15 simulation, produces very realistic results for many grain types. However^ the Cineon® 
application does not yield good performance for many high speed films because of the 
noticeable diagonal stripes the application produces for high grain size settings. Further, the 
Cineon® application fails to simulate grain with adequate fidelity when images become 
subject to prior processing, for example, such as when the images are copied or digitally 

20 processed. 

Another commercial product that simulates film grain is Grain Surgery™ from 
Visual Infinity Inc., which is used as a plug-in of Adobe ® After Effects ®. The Grain 
Surgery™ product appears to generate synthetic grain by filtering a set of random numbers. 
This approach suffers from disadvantage of a high computational complexity. 
15 Thus, a need exists for an efficient film grain simulation technique, which reduces the 

need for memory bandwidth, and computational effort, thus permitting film grain simulation 
in cost-sensitive high volume devices, such as set top boxes. 
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BRIEF SUMMARY OF THE INVENTION 

Briefly, in accordance with a preferred embodiment of the present principles, there is 
provided a method for simulating film grain in an image block of M x N pixels, where N and 
M are integers greater than zero. The method commences by first computing the average of 
the pixel values within the block of MxN pixels. A film grain block ofMxN pixels is 
selected from among a pool of previously established blocks containing film grain as a 
function of the average value of the image block and a random number. Each pixel in the 
selected film grain block is blended with a corresponding pixel in the image block. 

BRIEF DESCRIPTION OF THE DRAWINGS 

■ • 

FIGURE 1 depicts a block schematic drawing of an apparatus for generating pre- 
established film grain blocks for use in subsequent film grain simulation; and 

FIGURE 2 depicts a block schematic drawing of an apparatus in accordance with the 
present principles for simulating film grain on a pixel-by-pixel basis using the pre- 
established film grain blocks generated by the apparatus of FIG. 1. 
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DETAILED DESCRIPTION 
Introduction 

5 The method of the present principles simulates film grain in accordance with film 

grain information transmitted with an image to which simulated grain is blended. In 
practice, the transmitted image typically undergoes compression (encoding) prior to 
transmission via one of a variety of well-known compression schemes, such as the H.264 
compression scheme. With the transmitted image compressed using the H.264 compression 

10 scheme, transmission of the film grain information typically occurs via a Supplemental 

Enhancement Information (SET) message. Pursuant to contributions recently adopted by the 
standards body responsible for promulgating the H.264 standard, the SEI message can now 
include various parameters that specify different film grain attributes. 
Constraints on the film grain SEI message parameters 

■ 

15 

The method of the present principles imposes some constraints with regard to the number of 
parameters and their range of possible values allowed by the H.264 recommendation. 

_ * 

TABLE 1. provides a list of such parameters, including a description of their semantics and 
the constraints imposed by the present principles. 

>0 

TABLE I 



FILM GRAIN PARAMETER 


DESCRIPTION & CONSTRAINTS 


model_id 

* 


This parameter specifies the simulation 
model. It shall be 0, which identifies the 
film grain simulation model as frequency 
filtering. 


separate_colour_description_present_flag 


This parameter specifies if the color space in 
which the parameters are estimated is 
different from the color space in which the 
video sequence (where the film grain SEI 
message has been embedded) has been 
encoded. It shall be 0, which identifies the 
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• 


v^vyiwx apat^c iui nixn graJIl lilc- aalQw Ulan LUC 

encoded sequence. 


blehding^modejd 


This parameter identifies the blending mode 
used to blend the simulated film grain with 
the decoded images. It shall be 0, which 
correspond to an additive blending mode. 




iiiia paxaiiicicr luenuries me logaritximic 
scale factor used to represent the film grain 
parameters in the SEI message. It shall be in 
the range [0, 4] to ensure film grain 
simulation can be performed using 16-bit 

aTllillllCLllj. 


comp_modeI_present_flag[l] 


This parameter enables the transmission of 
film grain parameters for the Cb color 
component in the YCbCr color space. It 
shall be 0, since film grain simulation in 

tXuv/ina lo not 2>UppOIT€G. 


comp_model_present_flag[2] 


This parameter enables the transmission of 
film grain parameters for the Cr color 
component in the YCbCr color space. It 

CnQl 1 O CltiOA film frroin cimnliflrv-n ivi 

5iiaii uc u, amce rum groin simuiaiion in 
chroma is not supported. 


no Jntensity_intervals_minusl [0] 


This parameter defines the number of 
intensity intervals for which a specific set of 
parameters has been estimated. It shall be in 
the range [0, 7]. 


in tensityjnterval Jo werjbound [0] [i+1], 
intensity interval upper boundroilil 


These parameters define the boundaries of 

the luma intensitv intervals for which 

different film grain parameters are defined. 
The lower bound of interval i+1 must be 
greater than the upper bound of interval i 
because multigenerational film grain is not 
allowed. 
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num_jmodel_values_minusl[0] 


This parameter specifies the number of 

* 

model values present for each intensity 
interval in which the film grain has been 
modeled. It shall be in the range [0,4] 
because color correlation is not allowed. 


comp_model_value[0][i][0] 


This parameter represents the film erain 
intensity for each luminance intensity 
interval in which film grain has been 
modeled. It shall be in the range [0,255] to 
ensure film grain simulation can be 
performed using 16-bit arithmetic. 



In addition to the previous constraints, the present principles imposes that film grain SEI 
messages precede I pictures, and only one film grain SEI message can precede a particular I 
picture. (The presence in the bit stream of slice_type equal to 7 or nal_ref_idc equal to 5, 
5 indicates an I picture.) 



All the other parameters of the film grain SEI message have no constraint with respect to the 
standard specification. 

10 Bit-accurate implementation of film grain simulation 

Film grain simulation in accordance with the present principles occurs in a two-step 
process. First, generation of a pool of film grain blocks occurs during initialization, as 
described in greater detail with respect to FIG. 1. Thereafter, selected film grain portions are 

L5 added to each luminance pixel of each decoded picture as described with respect to FIG. 2. 

FIGURE 1 depicts an apparatus 10 in accordance with an illustrated embodiment of 
the present principles for generating a pool of film grain blocks for use in film grain 
simulation. The apparatus 10 typically generates a pool of 128 film grain blocks for each of 
as many as 8 different luminance intensity intervals. The SEI message field 

10 num ^intensity jntervals_minusl[0] indicates one less than the number of the luminance 
intensity intervals. 
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The apparatus 10 accomplishes film grain noise initialization using a specified 
uniform pseudo-random number polynomial generator 12 and using a specified list of 2048 
8-bit Gaussian distributed random numbers stored in a look-up table 14 The look-up table 
14 stores random numbers in 2's complement form in the range [-63, 63], The list of 
5 Gaussian random numbers appears in the Appendix. 

According to the bit-accurate specification of the present principles, generation of the 
film grain blocks begins with the lowest luminance intensity interval. The uniform random 
number generator 12 generates an index for the Gaussian random number list stored in the 
look-up table 14 using a primitive polynomial modulo 2 operator, x 18 + x? + x? + x 1 + 1. For 
10 ease of understanding, the term x(i, s) will indicate the i symbol of the sequence x, 

beginning with an initial seed s. The random number seed becomes reset to 1 upon the 
receipt of each film grain SEI message. 

To form an individual 8x8 film grain block, a random block generator 16 reads 8 
lines worth of 8 random numbers from the Gaussian random number look-up table 14. A 
15 random offset, from the random number generator 12, serves to access each line of 8 random 
numbers. Each line of the block produced by the block generator 16 is generated as 
following: 

index = x(i, 1) 

20 for n=0...7, B[i%8][n] = Gaussian Jist[(index + n)%2048] 
where i increments for each 8x1 block line. 

The 8x8 block of random values read by the generator 16 undergoes a transform, 
typically an integer Discrete Cosine Transform (DCT), performed by an Integer DCT 
25 transform block 18. After the DCT transform, the 8 x 8 random values undergo frequency 
filtering at a frequency filter 20 in accordance with the cut frequencies specified in the SEI 
message. Following frequency filtering, the 8x8 random values undergo an inverse DCT 
transform by an inverse integer DCT block 22. A first scaling block 24 scales the pixels on 
the top and bottom block lines as follows: 

30 

for n=0..7, B'[0][n] = (B[0] [n] + 1) » 1 
for n=0..7, B'[7][n] = (B[7] [n] + 1) » 1 
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This process continues until generation of a set of 128 film grain blocks for each luminance 
intensity interval. Following subsequent scaling by the second scaling block 26, the film 
grain blocks undergo storage in the film grain pool 28. 

5 Block and Pixel Operations to Simulate Film Grain 

FIGURE 2 illustrates an apparatus 200 in accordance with an illustrative embodiment 
of the present principles for simulating film grain on a pixel-by-pixel basis using the stored 
values in the film grain pool 28. The apparatus 200 includes a processing block 202 for 
10 creating an average of each 8x8 block of luma pixel values for comparison to the 
parameters intensity_intervaI_lower_bound[0][i] and 

ihtensity_intervaI_upper__bound[0][i] in the film grain SEI message to determine the 
correct luminance intensity interval for the current block. 

A selector block 204 selects a k* film grain block from the pool 28, using the.random 

L5 number generated by the uniform random number generator 16 from the polynomial modulo 
128 as the block index. Thus, the noise generator 16, which generates uniformly distributed 
random numbers using a polynomial for the initialization process described with respect to 
FIG. 1, finds application in the apparatus 200 of FIG. 2 to select film grain blocks, with the 
random number seed reset to 1 after the pool creation process. If the resulting block index is 

\0 identical to the previous one, the last bit of the index undergoes toggling. Such operation 
can occur using a bit-wise comparison and an XOR operator ( A ) as follows: 

previous_index = index 
index = x(k, 1) % 128 
:5 index A = (index = previous_index) 

■ 

Following block selection, a deblocking filter 206 deblocks the pixels on the right 
most column of the previously selected block and on the left most column of the current 
block. An adder 208 adds the deblocked film grain block to decoded luma pixels. (Since 
0 two horizontally adjacent blocks are required to perform deblocking, there is a 1 -block delay 
between the block selected in 204 and the block added in 208.) A clipper 210 clips the result 
within the range [0, 255] for display. Note that film grain noise addition only occurs to luma 
pixels. 
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Scaling of Cut Frequencies 

The parameters in the film grain SEI message of TABLE 1 assume the use of a 16 x 
5 16 DCT in the simulation process. In particular, horizontal and vertical high cut frequencies, 
provided by comp_model_value[0][i][l] and comp_model_value[0][i][2], and horizontal 
and vertical low cut frequencies, provided by comp_model_value[0][i][3] and 
comp_model_value[0][i][4], serve to filter the transform coefficients of a block of 16x16 
values. 

10 , 

In the illustrated embodiment, the use of 8 x 8 blocks will reduce complexity. 
Employing an 8 x 8 block transform using cut frequency parameters based on a 16 x 16 
transform implies that all the cut frequencies require scaling before the grain generation. 

The scaling of the cut frequencies occurs as follows: 

15 

comp_model_valueTO][i][j] = (comp_model_value [0][i][j] + 1) » 1 

where j is in the range [1,4]. Note that the scaling constitutes the equivalent of the integer 
division, rounded up to the nearest integer. 

20 

Integer Transform and Variance Scaling 

The transform used for the frequency filtering corresponds to an 8 x 8 integer 
approximation to the DCT, using the following transformation matrix: 

25 
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16-bit arithmetic can be used. The forward integer transformation of a given block of 
random noise is defined as: 

B = (((T 8 xB + 8) » 4)xT 8 r +8) » 4 
with 1 1 bits used for B . 



The inverse integer transform is defined as: 



L0 B = (T 8 r xBxT s +128)»8 
with 8 bits used for B. 

Following the inverse transform, the block B undergoes scaling as follows, assuming it is in 
the k luminance intensity interval, 

15 

val = B(i, j) * comp_model_value[0][k][0] 

B'(i, j) = (((val - (val »4) + 2 Iog2 -^ ale - factor " l ) » log2_scale__factor ) + 16) » 5 

where the operation (val - (val»4)) compensates the scaling of the integer transform; 
!0 log2_scale_factor, transmitted in the SEI message, scales comp_jmodel_value[0][k][0]; and 5 
scales the Gaussian numbers provided in the Appendix. 

Deblocking filter 206 

'5 As indicated, the film grain simulation apparatus of FIG. 2 includes a deblocking 

filter 206 for smoothing blocking artifacts resulting from the small size of the transform. In 
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the illustrated embodiment, the deblocking filter 206 takes the form of a 3-tap filter applied 
to all pixels bordering the 8x8 block left and right edges. Given a row of pixels belonging to 
two adjacent 8x8 blocks, the transition between blocks being located between pixels b and c, 



Block A 



Block B 



application of the deblocking occurs as follows: 



10 



b' = (a + (b « 1) + c) » 2 
c' = (b + (c « 1) + d) » 2 



15 



where b' and c 1 replace the value of the original pixels b and c, respectively. Deblocking of 
the left and right block edges occurs for every film grain block before addition to the 
decoded image. 

The foregoing describes a technique for simulating film grain in an image. 
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Appendix 

The list of the 2048 Gaussian distributed random numbers are: 
char Gaussian[2048] = { 

OxFB, 0x05, 0x33, OxFB, 0x14, OxEF, 0x06, OxlD, 0x26, 0x30, 0xD5, 0x01, 0x20, OxD9, 0x16, OxlB, . 
0xE7, OxOA 0x06, OxFB, 0xF6, OxF7, 0x10, OxCl, 0x08, OxFE, OxCC, 0x09, 0x09, 0x23, 0x17, OxFB, 
OxED, 0x15, OxFF, 0x25, OxDF, Oxl A, OxD3, 0x10, OxE9, OxOA, OxFF, 0xE5, 0x18, 0x00, OxE4, OxEC, 
0x00, Ox3C, OxCl, OxCB, OxE8, 0x04, 0x07, Ox3F, Ox3D, 0x36, 0x19, Ox3F, 0x00, 0x03, 0x38, 0x09, 
OxOE, 0x06, 0x26, 0x38, 0x28, 0x^2, OxCl, 0x37, OxE7, 0xF2, 0x01,.0xE8, OxF5, OxlD, OxF2, OxDC, 
0x05, 0x38, 0x21, 0x27, OxFF, 0xC7, 0xD5, OxFE, OxFE, 0x14, OxlD, 0xD8, 0x18, OxF3, OxFl, OxEF, 
OxCC, 0x19, 0x08, 0xF4, OxEF, OxFA, OxF9, OxCl, 0xE5, OxF5, OxE5, OxCl, OxC8, 0x02, 0xF4, OxDC, 
0x3F, Ox3F, OxFF, 0x14, 0x2B, OxEO, OxF9, OxlB, 0x09, 0x2D, OxD8, OxEO, OxEO, Oxl 1, OxFD, OxE5, 
0x31, OxFD, Ox2C, Ox3E, 0xF3, 0x2D, 0x00, Ox IF, OxlD, 0xF9, OxF5, 0x38, OxFO, Ox3A, 0x06, OxOC, 
0x19, OxF8, 0x35, OxFD, Oxl A 0x13, OxEF, 0x08, OxFD, 0x02, OxD3, 0x03, OxlF, OxlF, 0xF9, 0x13, 
OxEE, 0x09, OxlB, 0x08, 0xE7, 0x13, 0x10, OxEE, Ox3E, OxED, OxC5, 0x08, OxFl, 0x00, 0x09, 0x31, 
OxlE, 0x32, OxFA, OxDC, 0xF8 r OxE7, 0x31, 0x01, 0x01, OxlD, 0x10, OxFF, OxFF, 0x04, OxEC, OxCC, 
OxEE, 0x06, 0x3F, 0x07, OxCl, OxFl, OxD5, OxED, OxE5, 0x16, OxEC, 0x25, OxOB, OxF7, 0xF5, OxDD, 
0x25, 0xE6, 0x00, 0x10, OxEA, 0x08, 0xD2, OxlD, OxEO, OxDF, OxlB, OxCE, 0xF2, OxD5, OxEF, 0xD2, 
0x21, 0x02, OxDC, 0xE2, 0x2E, OxEB, 0x06, 0xF4; OxEE, OxCl, OxF8, 0x07, OxCl, OxlF, Oxll, OxOF, 
0x2E, 0x08, 0xE7, OxE3, 0x23, 0x26, 0x28, Ox3F, Ox3F, OxlE, 0x10, OxCC, 0xD2, 0x00, 0x00, 0x25, 
OxDE, 0x23, 0x3F, OxF7, OxC9, OxOE, OxOB, 0x07, 0x01, 0x13, Ox2D, 0x02, 0x14, 0x00, OxFE, 0x13, 
0x07, 0x38, 0xF2, OxEE, 0x19, 0x15, 0x35, OxOD, Ox3B, 0x03, OxD9, OxOC, OxDE, OxF6, 0x2E, OxFB, 
0x00, 0x09, 0x14, OxE7, 0x27, OxCl, OxEB, Ox3F, 0x08, 0x05, 0xF6, OxOF, 0xE7, OxOD, 0xD4, OxD3, 
OxED, OxF7, OxFC, OxOC, 0xC6, 0x23, 0xF4, OxEB, 0x00, 0x05, 0x2A OxCB, 0x13, OxFO, OxCl, 0x17, 
0x19, OxF4, OxF6, 0x16, 0x00, 0x07, OxEF, OxDE, 0x00, OxDC, OxOC, OxFD, 0x00, OxOE, OxFF, 0x16, 
0x10, OxFO, Ox3A, OxEA, 0x27, 0xF5, OxF8, OxCA, OxFB, OxDD, 0x2C, OxE9, OxOB, OxD3, 0x3B, OxEE, 
0x18, OxCl, OxlD, 0x10, OxD8, OxFB, OxF8, OxFD, 0x16, OxCl, 0xF9, 0x2C, Ox3F, 0x08, 0x31, OxED, 
OxFO, 0x12, 0x15, OxED, OxFl, OxF6, 0x34, OxF7, 0x09, 0x09, OxE3, OxFC, OxOF, 0x00, OxCl, 0x10, 
0x3F, OxD6, 0x25, OxOB, OxEC, OxE8, OxCl, OxCB, OxF9, 0x16, OxDB, 0x00, OxOE, OxF7, 0x14, OxDE, 
OxED, 0x06, Ox3F, OxFF, 0x02, OxOA, OxDC, 0xE3, OxCl, OxFF, OxFF, 0xE6, OxFE, OxC5, 0x2E, Ox3B, 
OxD8, OxE8, 0x00, 0x09, OxEA 0x21, 0x26, OxFA, OxF6, OxCl, Oxl 1, OxEC, OxlB, Ox3B, OxFE, 0xC7, 
OxF5, 0x22, 0xF9, OxD3, OxOC, OxD7, OxEB, OxCl, 0x35, OxF4, OxEE, 0x13, OxFD, OxFD, 0xD7, 0x02, 
OxD5, 0x15, OxEF, 0x04, OxCl, 0x13, 0x22, 0x18, OxEl, 0x24, OxE8, 0x36, 0xF3, 0xD4, 0xE9, OxED, 
0x16, Oxl 8, OxFF, OxlD, OxEC, 0x28, 0x04, OxCl, OxFC, 0xE4, OxE8, Ox3E, OxEO, 0x17, Oxl 1, Ox3A, 
0x07, OxFB, OxDO, 0x36, Ox2F, OxF8, OxE5, 0x22, 0x03, OxFA OxFE, 0x18, 0x12, OxEA, Ox3C, OxFl, 
OxDA 0x14, OxEA 0x02, 0x01, 0x22, 0x08, OxD9, 0x00, 0xD9, 0x02, 0x3F, 0x15, OxOD, Ox3F, OxCl, 
OxOD, 0xE5, OxF3, OxlB, 0x37, 0x17, 0x35, 0x00, OxDA, 0x00, OxlA, OxFC, 0xF5, OxEB, Ox3D, 0x36, 
Ox3F, 0x32, 0x21 , 0x17, 0x02, 0x00, 0x3D, OxFA, QxE5, OxFO, 0xE8, 0x2C, 0x20, OxCC, OxFE, Ox2F, 
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0xE6, OxlF, 0x16, OxOE, 0x17, 0x09, OxEF, 0x07, 0x14, 0x17, OxDO, OxF4, 0x2F, OxDB, Ox3F, 0xC7, 
Ox3F, OxDF, 0x00, OxF8, 0x19, OxDl, 0x17, 0x05, Oxll, OxEA, OxDB, 0x2C, OxCB, OxFC, 0xE4, 0xF2, 
OxCA, 0xF4, Ox3F, 0xE2, OxFA, 0x26, OxEA, 0x08, 0x09, 0x29, OxF5, 0x04, 0x3F, OxDF, OxlA, 0x01, 
OxOC, 0x06, 0x37, 0x15, 0xC8,.0xF5, 0x05, OxF4, 0x29, 0x21, OxFA, 0x25, OxC3, OxlD, Ox3F, OxFB, 
0x31, OxF7, OxlF, OxED, OxlA, 0x04, 0x03, OxlE, OxE5, 0x01, OxE4, 0x38, OxCC, OxE3, 0x01, OxFC, 
0xE9, 0x24, 0x2A, 0xE5, OxEF, 0x06, Ox3B, OxOD, 0x2E, OxDD, 0x06, OxCF, OxDD, OxF6, OxOE, 0x23, 
OxDl, 0x09, OxE6, 0x20, OxFA, OxEl, 0xF4, 0x20, 0x24, OxFC, 0x3F, 0x00, OxCl, 0x33, 0xF6, OxDC, 
0xC9, OxCD, OxFD, OxOE, OxEC, 0xF6, 0xE3, 0xF2, 0xF4, 0x09, OxFE, OxE7, 0x2F, OxE3, OxDl, OxEE, 
Oxl 1, 0x09, OxDE, Ox3F, OxF7, OxCl, 0xF5, 0xC5, OxE6, 0x12, 0x25, OxCl, 0x00, OxFB, 0xC5, 0xE6, 
OxF3, 0x13, 0x22, 0x08, 0x08, 0xC7, 0x2C, OxlF, OxOC, 0x12, 0xF5, 0x18, OxCE, OxFl, OxFC, OxDl, 
OxE6, 0x02, 0x2E, OxF5, OxE8, OxFC, 0x19, 0x01, OxDB, OxD4, OxFB, OxED, 0x3F, OxD5, 0xF5, 0x09, 
OxOA, 0x38, 0x25, 0x19, OxFl, 0x2E, OxEl, 0x03, OxFB, 0x17, 0x12, 0x32, OxEB, 0xF8, 0xE6, OxFD, 
OxEE, OxDA, OxFl, 0xF6, OxlF, OxOF, OxlF, OxOA, OxCl, OxOF, OxlF, 0x12, 0x33, 0xD6, OxFC, 0x26, 
0x27, OxlD, 0xD9, OxFD, Oxl 1, 0x04, 0x28, 0xF4, OxFC, 0x01, OxF8, 0x23, 0x3F, 0x29, OxD5, OxlB, 
0x09, 0xC5, OxC3, 0x12, 0x05, Ox3F, OxlC, OxE5, 0x38, 0x06, OxOC, 0x10, OxFA, 0xE9, OxOA, OxFA, 
0x02, OxlC, OxOD, OxOC, OxOC, OxFB, OxEE, 0x12, 0xD2, 0x26, 0x28, 0x04, 0x19, 0x06, 0x21, OxFA, 
0x00, 0x10, 0x16, OxDB, 0x10, OxED, 0xF5, OxE8, OxCl, OxF3, OxOF, OxFC, Oxll, 0x06, 0x23, 0x06, 
OxlC, 0x05, 0xE6, OxD6, OxlA, OxEA, OxEF, 0x00, Ox3F, 0x05, OxDF, OxEA, 0x17, 0xC7, 0x01, 0x05, 
OxlC, OxEF, Ox3B, OxF7, 0xE2, OxlA, OxE3, OxCl, OxE8, OxF5, 0x01, OxFE, 0x08, OxD8, OxFE, Ox3F, 
OxOC, 0x27, 0x21, OxlF, OxF4, 0x06, OxEO, OxEE, OxCl, 0xF2, OxOA, OxEl, 0x20, 0xE6, OxEC, 0x36, 
OxEl, 0x07, 0xF6, 0x06, OxOE, OxEl, OxOA, OxOD, 0x2F, OxEA, OxE3, 0xC6, OxFC, 0x27, OxE8, OxOB, 
OxEB, OxF8, 0x17, 0xE9, 0xC4, OxEF, 0xF2, OxE6, OxEA, OxOE, Ox3F, OxFA, 0x18, OxFC, OxCl, 0x25, 
OxF3, 0xF5, 0x2C, OxlD, 0x05, OxDl, 0x28, OxE3, OxlD, OxlE, 0xF4, 0x14, 0xD3, OxFF, 0xF6, OxE3, 
OxEA, OxE3, 0xF5,0xE6, 0x23,.0xF2, 0x21, OxFl, OxF5, 0x07, OxF8, OxDF, 0xF4, 0xF2, 0xE2, 0x17, 
0x12, 0x08, 0x07, OxEE, OxF5, OxFB, 0x04, OxF3, OxF7, OxlD, 0x16, 0xE8, 0xE9, OxFF, 0xF6, OxD8, 
OxOE, OxDF, OxCl, 0x25, 0x32, 0x02, OxF8, 0x30, Oxl 1, OxEO, 0x14, OxE7, 0x03, OxE3, OxOB, 0xE4, 
0xF7, 0xF4, 0xC5, OxDC, 0x2D, 0x07, 0xF9, 0x27, OxFO, 0xD9, OxCl, OxEF, 0x14, 0x26, OxD7, 0x00, 
OxlB, OxOB, OxDB, Ox3F, OxF8, 0xF6, 0x06, OxOF, OxlB, 0xC8, OxCl, 0x2C, OxlB, OxlE, 0x06, OxlB, 
OxFA, OxC8, 0xF9, OxOF, 0x18, OxDF, OxF8, 0x2D, OxFC, 0x00, OxOA, 0x22, OxDD, 0x31, OxF7, OxC8, 
0x20, OxD3, OxFC, OxFC, OxDD, Ox3F, 0x19, OxD8, 0xE8, OxOC, OxlE, 0xE2, OxC9, 0x03, OxEC, Ox3F, 
0x2B, OxEO, 0x35, OxCl, OxFE, Oxl 1, 0xF9, 0x14, OxE8, 0x06, 0x06, 0x24, OxCE, OxF3, 0x26, Ox3F, 
OxFD, OxCE, 0x2C, 0x12, Ox3C, 0x2C, 0xC2, OxE3, 0x06, 0xD2, 0xC7, OxOA, OxDF, 0xD5, OxDl, 0xC5, 
0x15, 0xF2, .OxFl, 0x08, 0x02, 0xE6, 0xE2, OxOA, OxEB, 0x05, OxDA, OxE3, 0x06, OxOE, 0x01, 0x03, 
OxDC, 0x13, OxE3, OxFB, 0x36, OxE6, 0x14, 0x21, OxFA, OxCl, OxCl, OxE8, OxOB, OxOE, 0x17, Oxl 1, 
0x2D, Oxll, OxFO, 0x39, 0xE7, OxFO, 0xE7, 0x2D, 0x03, OxD7, 0x24, 0xF4, OxCD, OxOC, OxFB, 0x26, 
Ox2A, 0x02, 0x21, OxD8, OxFA, OxF8, OxFO, OxE8, 0x09, 0x19, OxOC, 0x04, OxlF, OxCD, OxFA, 0x12, 
Ox3F, 0x38, 0x30, Oxl 1, 0x00, OxFO, OxE5, Ox3F, OxC3, OxFO, OxlE, OxFD, 0x3B, OxFO, OxCl, OxE6, 
OxEB, OxlF, 0x01, OxFE, OxF4, 0x23, 0xE4, OxFO, OxEB, OxEB, 0x10, 0xE4, OxCl, Ox3F, OxOC, OxEF, 
OxFB, 0x08, OxD8, OxOE, 0xE4, 0x14, OxCl, OxCl, OxOA, 0xE9, OxFB, OxEF, OxEl, 0xE7, OxFO, OxD8, 
0x27, OxDA, OxDC, 0x04, OxOD, OxDC, OxFC, OxDB, 0xD6, 0xD6, OxE4, OxOC, 0x27, OxFC, OxDO, Oxl 1, 
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OxEO, 0x04, 0xE3, 0x07, 0x00, OxEC, 0x10, 0xD5, OxiEA, 0x08, OxFF, OxFC, OxlD, 0x13, 0x05, OxCA, 
OxED, OxOB, 0x10, 0x08, 0xF2, 0x01, 0x19, OxCA, OxFE, 0x32, 0x00, 0x20, OxOB, 0x00, 0x3F, OxlE, 
0x16, OxOC, OxFl, 0x03, 0x04, OxFD, OxE8, 0x31, 0x08, 0x15, 0x00, OxEC, 0x10, OxED, 0xE6, 0x05, 
OxCA, OxF7, OxlC, OxCl, 0x22, OxOD, 0x19, 0x2E, 0x13, OxlE, 0xE7, 0x16, OxED, 0x06, 0x2A, Ox3C, 
OxOD, 0x21, 0x16, 0xC9, OxD7, OxFF, OxOF, 0x12, 0x09, OxEE, OxlD, 0x23, 0x13, OxDA, OxE9, OxlD, 
0xD9, 0x03, OxEl, OxEF, OxFA, OxlE, 0x14, OxCl, 0x23, OxFE, OxOB, 0xE5, 0x19, OxCl, 0x21, OxFE, 
OxEC, OxOE, OxEl, OxlD, OxFF, 0x00, OxF7, OxEA, 0xD2, OxD8, OxDO, 0xF9, 0xE6, OxFB, OxFB, OxDA, 
0x06, 0x00, 0x03, OxDF, OxCl, Ox3F, OxF3, OxOD, OxFA, 0x08, OxFA, OxF3, 0x00, 0x04, OxE9, OxFO, 
OxF9, OxOD, OxFl, OxE3, OxlD, 0x26, 0xC4, OxOD, 0x13, OxE5, OxEl, OxFl, 0xF6, OxEE, OxFl, OxFD, 
OxCl, OxF4, 0xE2, 0x23, OxCl, 0x38, OxCl, 0x3F, 0x2B, OxFD, 0x39, 0x36, OxlA, 0x2B, OxCl, 0x01, 
0x07, OxOB, 0x25, OxCC, OxE7, 0x01, 0x24, OxD8, 0xC9, OxDB, 0x20, 0x28, OxOC, OxlA, Ox3F, OxEA, 
0xE7, OxCD, OxEC, OxEO, 0xF2, 0x27, OxDF, 0x20, OxFO, OxFl, OxFD, Ox3F, 0x00, OxFA, 0xE7, 0x21, 
OxF9, 0x02, OxD2, OxOE, OxEF, OxFD, OxD3, OxE4, OxFF, 0x12, 0x15, 0x16, OxFl, OxDE, OxFD, 0x12, 
0x13, 0xE7, 0x15, OxD8, OxlD, 0x02, Ox3F, 0x06, OxlC, 0x21, 0x16, OxlD, OxEB, OxEB, 0x14, OxF9, 
OxC5, OxOC, 0x01, OxFB, 0x09, OxFA, 0x19, OxOE, 0x01, OxlB, OxE8, OxFB, 0x00, 0x01, 0x30, 0xF7, 
OxOE, 0x14, 0x06, 0x15, 0x27, OxEA, OxlB, OxCB, OxEB, 0xF7, Ox3F, 0x07, OxFB, 0xF7, 0xD8, 0x29, 
OxEE, 0x26, OxCA, 0x07, 0x20, 0xE8, 0x15, 0x05, 0x06, OxOD, OxOD, OxlE, OxlC, OxOF,.OxOD, 0x35, 
0xF7, OxlB, 0x06, 0x30, 0x02, OxFD, 0xE2, OxCD, 0x2F, 0x35, OxEB, OxlA, OxOD, 0xE9, OxFC, 0x34, 
0xE6, 0x17, 0x2C, 0x33, OxFO, 0x13, OxEF, OxlB, 0x19, 0x23, OxDl, OxEF, 0xD5, OxCB, 0xF7, OxFl, 
0x04, 0xF7, 0x27, 0xF9, 0x26, 0x02, 0XF7, OxCB, 0x2A, OxOA, OxEA, OxED, OxEC, 0x04, 0xF2, 0x25, 
0x17, OxDB, OxlE, OxCl, Ox3C, 0xC9, 0xE4, OxFl, 0x14, 0x03, 0x27, 0x25, 0x21, OxlC, 0x14, OxF4, 
OxOF, 0x12, 0xE9, OxEE, 0x15, OxDC, OxEE, OxlF, Ox3F, OxDE, OxE7, 0x2C, OxFO, OxE2, OxlD, 0xE5, 

• 

0x15, 0x07, 0x02, OxDF, 0x06, 0xD3, OxlF, OxOE, OxED, OxFF, 0x29, OxFF, OxED, 0xD6, OxD6, OxlC, 
Oxll, OxDE, 0xE2, OxOE, OxEE, OxDl, 0xD9, 0x02, OxOF, OxFE, OxFO, 0xD9, 0xF6, OxFC, OxDA, 0x16, 
0x03, 0xD2, OxDD, 0x20, 0x04, OxE8, Ox3F, OxDE, OxOC, OxFB, OxED, 0xC7, OxlF, OxCl, OxCE, 0x02, 
OxFl, 0x37, OxOB, OxE3, 0x20, OxCE, OxOD, OxEB, OxOA, OxE3, OxF3, OxDC, 0x01, 0xD2, 0x02, 0x3F, 
0x02, 0x25, OxD5, OxFC, OxEB, OxCE, Ox3F, 0x00, Ox3E, 0x2D, OxEl, 0x19, OxlC, 0x01, 0x28, OxCl, 
Ox3F, 0x27, Ox3F, 0xF2, OxOE, Ox3A, OxDB, OxF8, 0xE4, 0x34, 0x18, 0x16, OxOC, OxDD, 0x18, OxED, 
OxCB, OxOF, OxFO, 0x01, OxFB, 0x14, OxCl, 0x19, OxCC, OxEB, OxEE, 0x19, 0x00, 0x17, 0x2B, OxFC, 
0x26, OxOD, OxEC, 0xF4, 0x2D, 0x2B, 0xE5, 0x25, 0x05, 0x10, 0x26, OxlD, 0x3F, Ox3F, OxFD, OxDC, 
0x18, OxFO, OxCB, OxEF, 0x12, OxlC, OxlA, OxF8, OxFE, 0x29, OxlA, OxCB, OxlA, 0xC2, OxOE, OxOB, 
OxlB, OxEB, OxD5, OxF8, OxFD, 0x17, OxOB, OxFC, 0x00, OxFA, 0x37, 0x25, OxOD, 0xE6, OxEE, OxFO, 
0x13, OxOF, 0x21, 0x13, 0x13, OxEl, 0x12, 0x01, OxOA, OxFl, 0xE7, 0xF3, OxlA, OxED, 0xD5, OxOA, 
0x19, 0x39, 0x09, OxD8, OxDE, 0x00, 0xF9, OxE9, OxEA, OxFF, 0x3E, 0x08, OxFA, OxOB, 0xD7, 0xD7, 
OxDE, 0xF7, OxEO, OxCl, 0x04, 0x28, OxE8, OxlE, 0x03, OxEE, OxEA, OxEB, OxlC, 0xF3, 0x17,0x09, 
OxD6, 0x17, OxFA, 0x14, OxEE, OxDB, 0xE2, 0x2A, 0xD9, OxCl, 0x05, 0x19, 0x00, OxFF, 0x06, 0x17, 
0x02, 0x09, 0xD9, OxE5, OxF3, 0x20, OxDD, 0x05, OxCB, 0x09, OxF8, 0x05, OxFl, OxlF, OxE5, 0x12, 
0x25, OxF8, Ox3F, OxDC, OxFO, 0xF2, 0xC5, 0x34, 0x21, 0x35, OxCD, OxCC, 0x23, OxlE, 0x01, OxOB, 
OxFF, 0x10, OxFE, OxF9, OxDF, 0xF9, 0xF5, OxE5, 0x07, OxEl, 0x25, OxlC, OxC9, 0x00, 0x29, OxF3, 
OxOA, 0x25, OxED, OxF8, OxFB, 0x20, OxF8, OxCl, OxE5, OxEO, OxOF, 0x2F, Ox3A, 0x01, OxC8, OxFD, . 
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OxCA, OxEl, 0x30, 0x04, 0x19, 0x03, 0x25, 0xF3, 0x24, 0x38, OxEE, 0xC9, 0x2F, 0xE7, OxOB, OxFA, 
OxF7, OxlB, OxOA, OxOB, 0x2D, 0x2D, OxOB, OxE8, 0x08, OxDB, OxOB, 0x04, 0xE8, OxDO, OxEE, 0x18, 
OxEF, Oxll, OxCl, 0xD6, 0x15, Ox3F, 0xF5, 0xF4, 0x2A, 0x29, OxEF, OxFO, OxFA, 0x36, 0x33, OxED, 
0x19, OxDF, Oxll, 0x09, 0xF5, 0x18, OxFl, 0x3F, 0x14, OxOC, 0xD2, OxFF, OxFF, 0x34, 0x01, 0xE4, 
OxF8, 0x03, Ox3F, 0xF8, 0x3E, 0x21, 0x22, 0xE2, OxOF, OxEF, Oxl A, 0xE4, OxF5, 0x08, 0x15, OxEF, 
0xF3, 0xE4, OxDF, 0xF6, OxFC, 0xE8, 0x21, 0x06, 0x20, 0x02, 0x17, OxlB, Ox3F, OxDB, 0x16, 0x2C, 
OxEO, OxFA, OxDA, OxD8, 0xD3, OxOB, OxOE, 0x10, OxED, 0xD5, OxFO, 0x30, 0xD3, 0x13, 0x04, OxEl, 
OxFF, OxFB, Ox3F, 0xE8, OxEE, 0xE5, OxOB, OxEF, OxEF, OxE6, 0x2C, OxD3, 0x00, 0x18, 0x26, OxFE, 
OxCl, 0x08, 0x16, OxEE, OxDC, 0x00, 0xE4, OxF7, OxDC, OxOE, 0x2E, OxlD, 0x18, OxOA, 0x08, 0x37, 
0xC9, 0x10, 0xD7, 0x17, 0x17, OxFB, Oxl 1, OxD5, 0x15, OxlC, OxDO, Ox3F, 0xF8, 0x00, 0x00, OxED, 
OxCl, OxFF, 0x00, OxlF, 0x2E, 0x00, 0x12, OxEO, 0xE2, OxF7, 0x13, OxCl, OxlC, 0x18, OxF8, Ox3F, 
0x2C, OxEB, OxCA, OxE7, OxF8, 0x03, OxEE, 0x22, 0x17, OxF9, 0x35, 0x14, OxlC, 0x03, 0x09, 0x03, 
0x01, Ox2B, 0xD4, 0xD2, OxF8, 0xF6, 0xF5, 0x06, 0x03, OxFE, OxDA, 0xD3, OxFF, 0x03, OxEF, OxFE, 
0x09, 0x01, 0xC9, 0x02, OxDF, 0xD8, 0x3C, OxF7, OxFO, OxEE, 0xD6, Ox3F, 0x21, 0x16, 0x08, 0x17 

}; 
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CLAIMS 



1 1 . A method for simulating film grain in an image block of M x N pixels, where 

2 N and M are integers greater than zero, comprising the steps of: 

3 computing the average of the pixel values within the block of M x N pixels; 

4 selecting a film grain block of M x N pixels from among a pool of previously 

5 established blocks containing film grain as a function of the average value of the image 

6 block and a random number; and 

7 blending each pixel in the selected film grain block with a corresponding pixel in the 

8 image block. 

1 2. The method according to claim 1 further including the step of accessing a 

2 look up table containing random numbers to obtain the random number. 

1 3. The method according to claim 2 further comprising the step of populating 

2 the look-up table in advance of film grain simulation with random numbers generated by a 

3 random number generator. 



4. A method for creating a block of M x N pixels with film grain, where N and 
M are integers greater than zero, comprising the steps of: 

receiving film grain information that includes at least one parameter that specifies an 
attribute of the film grain to appear in the block; 

creating a block of M x N random values selected from a previously established list 
of Gaussian random numbers; 

computing an Discrete Cosine Transform of the M x N block of random numbers; 

filtering the M x N coefficients resulting from the Discrete Cosine Transform by at 
least one parameter in the received film grain information; 

computing an Inverse Discrete Cosine Transform of the filtered set of coefficients; 

scaling all the pixel values in the block as indicated by one parameter in the received 
film grain information; and 

storing the created block of film grain into a pool of film grain blocks. 
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1 5. The method according to claim 4 further comprising the step of performing an 

2 integer approximation of a Discrete Cosine Transform (DCT) and the Inverse Discrete 

3 Cosine Transform (IDCT) to reduce complexity. 

1 6. The method according to claim 4 further comprising the step of scaling top . 

2 . and bottom edges of the created film grain block to hide block edges. 

1 7. The method according to claim 4 wherein the step of receiving the film grain 

2 information further comprises the step of decoding a Supplemental Enhancement 

3 Information message containing the at least one parameter. 

1 8. Apparatus for simulating film grain in an image block of M x N pixels, where 

2 N and M are integers greater than zero, comprising: 

3 means for computing the average of the pixel values within the block of M x N 

4 pixels; 

5 means for selecting a film grain block of M x N pixels from among a pool of 

6 previously established blocks containing film grain as a function of the average value of the 

7 image block and a random number; and 

8 means for blending each pixel in the selected film grain block with a corresponding 

9 pixel in the image block. 

1 9. . The apparatus according to claim 8 further a look up table containing random 

2 numbers to obtain the random number. 

1 10. The method according to claim 9 where the look-up table is populated in 

2 advance of film grain simulation with random numbers generated by a random number 

3 generator. 



11. An apparatus for creating a block of M x N pixels with film grain, where N 
and M are integers greater than zero, comprising: 

means for receiving film grain information that includes at least one parameter that 
specifies an attribute of the film grain to appear in the block; 
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means for creating a block of M x N random values selected from a previously 
established list of Gaussian random numbers; 

means for computing an Discrete Cosine Transform of the M x N block of random 
numbers; 

means for filtering the M x N coefficients resulting from the Discrete Cosine 
Transform by at least one parameter in the received film grain information; 

means for computing an Inverse Discrete Cosine Transform of the filtered set of 
coefficients; 

means for scaling all the pixel values in the block as indicated by one parameter in 
the received film grain information; and 

means for storing the created block of film grain into a pool of film grain blocks. 

12. The apparatus according to claim 1 1 further comprising means for performing 
an integer approximation of a Discrete Cosine Transform (DCT) and the Inverse Discrete 
Cosine Transform (IDCT) to reduce complexity. 

13. The apparatus according to claim 11 further comprising the means for scaling 
top and bottom edges of the created film grain block to hide block edges. 

14. The apparatus according to claim 11 wherein means for receiving the film 
grain information further comprises means for decoding a Supplemental Enhancement 
Information message containing the at least one parameter. 
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